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Abstract 

(f) ' In this article, we study the vector and axialvector B c mesons using the QCD sum rules, and 

make reasonable predictions for the masses and decay constants, then calculate the leptonic 
decay widths, the present predictions can be confronted with the experimental data in the 
future. 






X 



PACS number: 14.40.Pq, 12.38.Lg 
Key words: B c mesons, QCD sum rules 

1 Introduction 



i-^ ' In 1998, the CDF collaboration observed the pseudoscalar bottom-charm B c mesons through the 

decay modes B^r —> J/ipi^X and B^r — > J/tpi^Di in pp collisions at the energy \/s = 1.8 TeV at 
Q_|i the Fermilab Tevatron, the measured mass is M Bc = (6.40 ± 0.39 ± 0.13) GeV [L. In 2007, the 

CDF collaboration observed the pseudoscalar B c mesons with a significance exceeds 8 a through 
the decay modes B^r — > J/tpir^ in pp collisions at the energy y/s — 1.96 TeV using the Collider 
Detector at Fermilab (CDF II), the measured mass is M Bc = (6275.6 ±2.9 ±2.5) MeV [2J. In 2008, 
the DO collaboration reconstructed the B^r — > J/ipn^ decays and observed the pseudoscalar B c 
mesons with a significance more than 5 a, the measured mass is Mb c — (6300 ± 14 ± 5) MeV [3]. 

I/-n ■ Now the average value listed in the Review of Particle Physics is Mb c = (6.277 ± 0.006) GeV [4]. 

Other B c mesons, such as the scalar, vector, axialvector, tensor B c mesons, have not been observed 

^sO ■ yet, but they are expected to be produced at the Large Hadron Collider (LHC) in the future [S]. 

rf\ ' The heavy quarkonium states and triply-heavy baryon states play an important role both in 

studying the interplays between the perturbative and nonperturbative QCD and in understanding 
the heavy quark dynamics due to the absence of the light quark contaminations. The bottom- 
charm quarkonium states B c , which consist of the heavy quarks with different flavors, are of spe- 
cial interesting. The ground states and the excited states lying below the BD, BD* , B*D, B*D* 
thresholds cannot annihilate into gluons, and therefore are more stable than the corresponding 
charmonium and bottomonium states, and would have widths less than a hundred KeV j6j- The 
excited states can undergo radiative or hadronic transitions to the ground state pseudoscalar B c 
mesons, which decay weakly. There have been several theoretical works on the mass spectroscopy 
of the B c mesons, such as the relativized (or relativistic) quark model with an special phenomeno- 
logical potential [5] [7J [5] [2], the nonrelativistic quark model with an special phenomenological 
potential [TO] [TTJ [12] , the semi-relativistic quark model using the shifted large- N expansion [13] , 
the perturbative QCD [14], the nonrelativistic renormalization group [15], the lattice QCD [161 ITT] . 
etc. 

The QCD sum rules is a powerful theoretical tool in studying the heavy quarkonium states 
[T51Q33, an d t ne existing works focus on the S'-wave heavy quarkonium states J/ip, rj c , T, 77^, and 
the P-wave spin-triplet heavy quarkonium states \cj, Xbj, j = 0, 1, 2 [19, 20]. The pseudoscalar B c 
mesons have been studied by the full QCD sum rules [21] [22] [23] [24] and the potential approach 
combined with the QCD sum rules [TT] [25] [26] , while the vector B c mesons (B*) have been studied 
by the full QCD sum rules [23] [24], and the axialvector B c mesons have not been studied yet. 
In Ref.[23], Colangelo, Nardulli and Paver took the leading-order 0(1) approximation, obtained 
the values Mb c ~ 6.35 GeV and fs- ~ Jb c = (360 ± 60) MeV, and did not present the value 
Mb«. In Ref . [24] . Narison took into account the next-to- leading-order 0(a s ) contributions by 
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assuming that one quark had zero mass, obtained the values Mb* — Mb* = (1.53 ± 0.18) GeV, 

J B * = 4^F and 7s; = 14 - 0± L0 ' the mass M b* = m b* + (1.53 ±0.18) GeV is much larger than 
other theoretical calculations [HI El [HI El OH [HI [HI EH] • Those studies based on the QCD sum rules 
were preformed before the pseudoscalar B c mesons were observed by the CDF collaboration, the 
predictions should be updated. Now we can take the experimental data as guides to choose the 
suitable Borel parameters and continuum threshold parameters. Naively, we expect that the masses 
of the pseudoscalar, vector and axialvector B c mesons have the hierarchy: M B r -\ < M Bc n-\ < 
Mb c (i+), the 0~, l - and 1 + denote the spin-parity J p . Furthermore, the calculations based on 
the nonrelativistic renormalization group indicate that M Bc (\-) — M Bc (q-) = (50 ± 17^^) MeV 
[T5] . In this article, we recalculate the masses and decay constants of the vector and axialvector 
B c mesons with the QCD sum rules by including the next-to- leading-order 0(a s ) contributions, 
and make predictions for the leptonic decay widths. The decay constants are basic parameters in 
studying the exclusive processes of the B c mesons. 

The article is arranged as follows: we derive the QCD sum rules for the masses and decay 
constants of the vector and axialvector B c mesons in Sect. 2; in Sect. 3, we present the numerical 
results and discussions; and Sect. 4 is reserved for our conclusions. 

2 QCD sum rules for the vector and axialvector B c mesons 

In the following, we write down the two-point correlation functions IT^,,^) in the QCD sum rules, 

IW(p) = i|crW^<0|T{^(a04(0)}|0>, (1) 

Jj(x) = c(x)^b{x) , 

tf(x) = cixfrwKx) , (2) 

where J^{x) = J^(x), j£(x), the vector and axialvector currents JY (x) and J^(x) interpolate the 
vector and axialvector B c mesons, respectively. 

We can insert a complete set of intermediate hadronic states with the same quantum numbers 
as the current operators J^{x) into the correlation functions H fll ,(p) to obtain the hadronic repre- 
sentation [TBI [H] . After isolating the ground state contributions from the vector and axialvector 
B c mesons, we get the following result, 

f\J A B c ( „ , PiiPv 

-2 
B c 

= n ( p)(-^ + M^+..., (3) 

where the decay constants fg c are defined by 

(0\UQ)\B c (p)) = f Bc M Bc e,, (4) 

and the e^ are the polarization vectors of the vector and axialvector B c mesons. The Tl(p) can be 
expressed in the following form through the dispersion relation, 

r°° fl Ml 5(s - Ml ) 

n(p) = / ds JB < Bc _ 2 — ^ + • • • . (5) 

J (m b +m c ) 2 s P 



M L ~ P \ P 



Now, we briefly outline the operator product expansion for the correlation functions 11^,, (p) 

l V/ 



in perturbative QCD. We contract the quark fields in the correlation functions n M /, (p) (here we 






Figure 1: The typical diagrams we calculate in the operator product expansion. 



add the indexes V and A to denote the vector and axialvector currents respectively) with Wick 
theorem firstly, 



KM 



n£,(p) 
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d A xe^ x Tr { ltllb B l0 {x) lvlb C 3l (-x)} , 



where the Bij(x) and Cij(x) are the full heavy quark propagators, and can be written as Sij(x) 
collectively, 
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and t n — 4-, the A™ is the Gell-Mann matrix, the i, j are color indexes, D a = d a — ig s G™t n , and 
the (g 2 GG) = (fi^G^G™^) is the gluon condensate [19]; then complete the integrals both in the 
coordinate space and in the momentum space; finally obtain the correlation functions Yi^ v {p) (or 
n(p)) at the level of the quark-gluon degrees of freedom. Furthermore, we calculate the 0(a s ) 
corrections to the perturbative terms originally. The typical Feynman diagrams we calculate in 
the operator product expansion are shown in Fig.l. 

Once analytical expressions of the correlation functions at the quark level are obtained and 



expressed as /, 



ds- 



P±(s) + ^p¥(s) 



through the dispersion relation, then we can 



take the quark-hadron duality and perform the Borel transforms with respect to the variable 
P 2 = —p 2 to obtain the QCD sum rules, 
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X(a, b,c) — a 2 + b 2 + c 2 — 2ab — 2ac — 26c, the p" (s) and p_ * (s) are the QCD spectral densities 
of the vector and axialvector B c mesons respectively, the T is the Borel parameter, and the So 
are the continuum threshold parameters. In this article, we take into account the contributions 
from the gluon condensate ( asGG ) and the four quark condensates (q~7 f j,t a qqj ,1 t a q), and neglect 
the contributions from the three-gluon condensate due to its small value. In calculations, we have 
used the equation of motion, D v Ct v — ^2„— u ^ s 9sq"f^t a q, and assumed vacuum saturation for the 
four quark condensates and taken the approximation (ss) = (qq) . 

We calculate the imaginary parts ° p^'(s) of the 0(a s ) corrections with the Cutkosky's 
rule firstly [37] , then use the dispersion relation to obtain the full 0(a s ) corrections. The calcu- 
lations should be carried out in the £)-dimension, because there exist divergences in the integrals 
ds^3-Tp^ s (s). In calculations, we observe that the divergences (endpoint diver- 



»00 r°° 
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gences) are canceled out with each other in the imaginary parts °' p± s (s), the imaginary parts 
° s w p^. 3 (s) have no divergences themselves, divergences only occur in the integrals over s. The 
integrals as ^ J> f. " 2 dsp± s (s) exp (— ^J are convergent and we can take the dimension D = A 

2 2 

safely The terms like log ^ and log 7 -^*t do not appear. In this article, we carry out the calculations 
by assuming m c ^ 0, the perturbative 0(a s ) corrections to the corresponding correlation functions 



II(p) of the heavy-light currents cannot be reduced by simply setting m c = [28], as there appear 
additional endpoint divergences, we should take care of the subtleties in performing the dispersion 
relation and recalculate the Feynman diagrams [27) . There exist three typical energy-scales m 2 , 
m 2 and p 2 , it is difficult to use the usual procedure of the dimensional regularization (through 
Feynman parameters) to carry out the two-loop divergent integrals. 

We can eliminate the decay constants /b c (i=f) & n d obtain the QCD sum rules for the masses of 
the vector and axial vector B r mesons, 
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then use the resulting masses as input parameters to obtain the decay constants /b c (it) from 
Eq.(7). 

3 Numerical results and discussions 

The mass of the pseudoscalar B c meson is Mb c = (6.277 ± 0.006) GeV from the Particle Data 
Group [4] , while the calculations based on the nonrelativistic renormalization group indicate that 
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parameters as s 
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(50 ± 17±i|) MeV 

= (41 - 47) GeV : 



[15. We can tentatively take the continuum threshold 
(46 — 54) GeV , and search for the ideal 



anci So / 1 -\ 



values, where we have assumed that an additional P-wave results in mass-shift about 0.5 GeV and 
the energy gap between the ground states and the first radial excited states is about 0.5 GeV. 

The vacuum condensates are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 
at the energy scale \x = 1 GeV [29] . The quark condensate evolves with the renormalization 



group equation, (qq)(fi 2 ) = (qq)(Q 2 ) 
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The value of the gluon condensate 



ia 3 GG\ 



has 



been updated from time to time, and changes greatly [20) . we use the recently updated value 
(s«|£) = (0.022 ±0.004) GeV 4 [3D]. 

In this article, we take the MS masses m c (ml) = (1.275 ± 0.025) GeV and m b {ml) = (4.18 ± 
0.03) GeV from the Particle Data Group [4|. Furthermore, we take into account the energy-scale 
dependence of the MS masses from the renormalization group equation, 
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or cb system is fi = 



'Bc(O-) 



= 5, 4 and 3, respectively [3]. The typical energy scale in the bc 
(m b + m c ) 2 rj (2.6 ± 0.1) GeV, where the spin-parity 0~ denotes 



the observed pseudoscalar meson B c [3], we can take n/ = 4. In Fig. 2, we plot the threshold 
(m b + m c ) 2 with variations of the energy scales. From the figure, we can see that the threshold 
(■nib + m c ) 2 decreases quickly with increase of the energy scale, the energy scale should be larger 
than 1.7 GeV for the be or cb system. The typical energy scale \i = (2.5 — 2.7) GeV leads to a 
rather large integral range so — (m& + m c ) 2 , it is a reasonable choice. 

In calculations, we observe that the ideal parameters are T 2 = (4.2—6.2) GeV 2 [(4.6—6.6) GeV 2 ] 
and sq = (47 ± 1) GeV 2 [(54 ± 1) GeV 2 ] for the vector [axialvector] B c mesons, the corresponding 
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Figure 2: The energy scale dependence of the threshold (m& + m c ) 2 , where 39.4 GeV 2 is the 
squared mass of the pseudoscalar meson B c . 





T 2 (GcV 2 ) 


5 (GcV z ) 


pole 


M(GeV) 


/s.(GcV) 


B c (l-) 


4.2-6.2 


47 ± 1 


(51 - 88)% 


6.34 ±0.09 


0.79 ±0.10 


Sc(l + ) 


4.6-6.6 


54 ± 1 


(51 - 85)% 


6.74 ±0.10 


0.81 ±0.11 



Table 1: The Borel parameters, continuum threshold parameters, pole contributions, masses and 
decay constants of the vector and axialvector B c mesons. 



pole contributions and the resulting masses and decay constants are presented in Table 1 and 
Fig. 3. From the table, we can see that the energy scale and the Borel parameters are of the 
same order, /i 2 = 0(T 2 ). The threshold parameters and the predicted masses satisfy the relations 

0.5 GeV and 



M F 



M F 



,, ! , x,j.b c {\~) ~ u.oijev aim .i a B (] — iub c (i+) ~ 0.6 GeV, which are compatible with 

our naive expectation that the energy gap between the ground state and first radial excited is 
about (0.5 — 0.7) GeV. From the table, we can also see that the uncertainties of the masses are 
much smaller than that of the decay constants, because the masses are obtained from a ratio, see 
Eq.(10), some uncertainties in the numerator and denominator are canceled out with each other. 
The calculations based on the nonrelativistic renormalization group indicate that Mg/i-\ — 

— 63 MeV based on 
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Bc(O-) 



(50 ± 17+ 1 2) MeV H3> the Present prediction M Ba 



(i- 



M 
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the MS masses plus 0(a s ) corrections is satisfactory, here we take the central value. 

In Table 2, we present the theoretical values of the vector and axialvector B c mesons from the 
relativized (or relativistic) quark model with an special potential (6] El HI IS] , the nonrelativistic 
quark model with an special potential [El HU H2] , and the lattice QCD [16]. From the Table, we can 
see that the present predictions are consistent with those values. In Table 3, we present the values 
of the decay constants of the vector and axialvector B c mesons from the relativistic quark model 
with an special potential [7] , the nonrelativistic quark model with an special potential [101 111! 112] , 
the light- front quark model [3TJ [32] , the Bethe-Salpeter equation [33] , and field correlator method 
[34] . Compared with those theoretical calculations f Bc (i-) = (380 - 520) MeV [7] IT0 | HT ] [T2 ] 
l3T| [32] 133] [34], the present prediction /s c (i-) = (790 ± 100) MeV is somewhat larger, while the 
present prediction /b c (i+) = (810 ± 110) MeV is much larger than the value 160 MeV from the 
Bethe-Salpeter equation [33) . At present time, it is difficult to say which value is superior to 
others. 
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Table 2: The masses of the vector and axialvector B c mesons from different theoretical approaches, 
the unit is GeV. 
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Table 3: The decay constants of the vector and axialvector B c mesons from different theoretical 
approaches, the unit is MeV. 

The leptonic decay widths Ti Pe of the vector and axialvector B c mesons can be written as, 
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where £ = e,n,T, the Fermi constant Gf = 1.16637 x 10 5 GeV 2 , the CKM matrix element 



V cb = 40.6 x 10~ 3 , the masses of the leptons m e = 0.511 x 10~ 3 GeV, m M = 105.658 x 10~ 3 GeV 
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m T = 1776.82 x 10~ 3 GeV 'A',. We can use the masses and decay constants of the vector and 
axialvector B c mesons in Table 1 to obtain the leptonic decay widths, 
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where the uncertainties originate from the uncertainties of the masses and decay constants, respec- 
tively. Such tiny leptonic decay widths maybe escape experimental detections. 

4 Conclusion 

In this article, we study the vector and axialvector B c mesons with the QCD sum rules, and make 
predictions for the masses and decay constants, which can be taken as basic input parameters in 
studying revelent physical physics, such as the semi-leptonic and non-leptonic decays, then calculate 
the leptonic decay widths, the present predictions can be confronted with the experimental data 
in the future at the LHC. 
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Figure 3: The masses and decay constants of the vector and axialvector B c mesons with variations 
of the Borel parameters T 2 , where (I) and (II) denote the vector and axialvector B c mesons, 
respectively. 



References 

[1] F. Abe et al, Phys. Rev. D58 (1998) 112004; F. Abe et al, Phys. Rev. Lett. 81 (1998) 2432. 

[2] T. Aaltonen et al, Phys. Rev. Lett. 100 (2008) 182002. 

[3] V. M. Abazov et al, Phys. Rev. Lett. 101 (2008) 012001. 

[4] K. Nakamura et al, J. Phys. G37 (2010) 075021. 

[5] C. H. Chang, Y. Q. Chen, G. P. Han and H. T. Jiang, Phys. Lett. B364 (1995) 78; K. 
Kolodziej, A. Lcikc and R. Ruckl, Phys. Lett. B355 (1995) 337; K. Cheung and T. C. Yuan, 
Phys. Rev. D53 (1996) 1232; K. Cheung and T. C. Yuan, Phys. Rev. D53 (1996) 3591; I. 
P. Gouz, V. V. Kiselev, A. K. Likhodcd, V. I. Romanovsky and O. P. Yushchenko, Phys. 
Atom. Nucl. 67 (2004) 1559; C. H. Chang and X. G. Wu, Eur. Phys. J. C38 (2004) 267; A. 
V. Berczhnoy, A. K. Likhoded and A. A. Martynov, Phys. Rev. D83 (2011) 094012. 

[6] S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189; S. Godfrey, Phys. Rev. D70 (2004) 
054017. 



[7; 

[9 
[10 
[11 

[12 
[13 

[14 
[15 

[16 

[17; 

[18 

[19 
[20 

[21 

[22 
[23 
[24" 



D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D67 (2003) 014027. 
S. N. Gupta and J. M. Johnson, Phys. Rev. D53 (1996) 312. 

J. Zcng, J. W. Van Ordcn and W. Roberts, Phys. Rev. D52 (1995) 5229. 

L. P. Fulchcr, Phys. Rev. D60 (1999) 074006. 

S. S. Gershtcin, V. V. Kiselev, A. K. Likhodcd and A. V. Tkabladze, Phys. Rev. D51 (1995) 
3613; S. S. Gershtein, V. V. Kiselev, A. K. Likhoded and A. V. Tkabladze, Phys. Usp. 38 
(1995) 1. 

E. J. Eichten and C. Quigg, Phys. Rev. D49 (1994) 5845. 

S. M. Ikhdair and R. Sever, Int. J. Mod. Phys. A19 (2004) 1771; S. M. Ikhdair and R. Sever, 
Int. J. Mod. Phys. A20 (2005) 6509; S. M. Ikhdair and R. Sever, Int. J. Mod. Phys. A20 
(2005) 403. 

N. Brambilla and A. Vairo, Phys. Rev. D62 (2000) 094019. 

A. A. Pcnin, A. Pineda, V. A. Smirnov and M. Steinhauser, Phys. Lett. B593 (2004) 124. 

C. T. H. Davies et al, Phys. Lett. B382 (1996) 131. 

E. B. Gregory et al, Phys. Rev. Lett. 104 (2010) 022001. 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385. 

L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1. 

S. Narison, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 17 (2002) 1. 

E. Bagan, H. G. Dosch, P. Gosdzinsky, S. Narison and J. M. Richard, Z. Phys. C64 (1994) 
57. 

M. Chabab, Phys. Lett. B325 (1994) 205. 

P. Colangclo, G. Nardulli and N. Paver, Z. Phys. C57 (1993) 43. 

S. Narison, Phys. Lett. B210 (1988) 238. 



10 



[25] V. V. Kiselev, A. K. Likhoded and A. I. Onishchenko, Nucl. Phys. B569 (2000) 473. 

[26] V. V. Kiselev and A. V. Tkabladze, Phys. Rev. D48 (1993) 5208. 

[27] S. Bauberger, M. Bohm, G. Weiglein, F. A. Berends and M. Buza, Nucl. Phys. Proc. Suppl. 
37B (1994) 95; S. Bauberger, F. A. Berends, M. Bohm and M. Buza, Nucl. Phys. B434 
(1995) 383. 

[28] A. Khodjamirian and R. Ruckl, Adv. Ser. Direct. High Energy Phys. 15 (1998) 345. 



[29] P. Colangelo and A. Khodjamirian, |hep-ph/0010175 ; B. L. Ioffe, Prog. Part. Nucl. Phys. 56 
(2006) 232. 

[30] S. Narison, Phys. Lett. B693 (2010) 559; S. Narison, Phys. Lett. B706 (2012) 412; S. Narison, 
Phys. Lett. B707 (2012) 259. 

[31] H. M. Choi and C. R. Ji, Phys. Rev. D80 (2009) 054016. 

[32] C. W. Hwang, Phys. Rev. D81 (2010) 114024. 

[33] G. L. Wang, Phys. Lett. B650 (2007) 15; G. L. Wang, Phys. Lett. B633 (2006) 492. 

[34] A. M. Badalian, B. L. G. Bakker and Yu. A. Simonov, Phys. Rev. D75 (2007) 116001. 



11 



